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ABSTRACT 

The measured redshift (z) of an astronomical object is a combination of Hubble recession, 
gravitational redshift and peculiar velocity. The line of sight distance to a galaxy inferred 
from redshift is affected by the peculiar velocity component of galaxy redshift, which is ob¬ 
served as an anisotropy in the correlation function. This anisotropy allows us to measure the 
linear growth rate of matter (fas). We measure the fag at 2 ; = 0.57 using the CM ASS sample 
from Data Release 11 of Sloan Digital Sky Survey III (SDSS III) Baryon Oscillations Spec¬ 
troscopic Survey (BOSS). The galaxy sample consists of 690,826 massive galaxies in the 
redshift range 0.43-0.7 covering 8498 deg^. Here we report the first simultaneous measure¬ 
ment of fag and background cosmological parameters using Convolution Lagrangian Pertur¬ 
bation Theory (CLPT) with Gaussian streaming model (GSRSD). We arrive at a constraint of 
fag = 0.462 ± 0.041 (9% accuracy) at effective redshift (z = 0.57) when we include Planck 
CMB likelihood while marginalizing over all other cosmological parameters. We also measure 
bag = 1.19±0.03, H{z = 0.57) = 89.2±3.6kms-i Mpc-^ andDA(^ = 0.57) = 1401±23 
Mpc. Our analysis also improves the constraint on = 0.1196 ± 0.0009 by a factor of 
3 when compared to the Planck only measurement(r2c/i^ = 0.1196 ± 0.0031). Our results 
are consistent with Planck ACDM-GR prediction and all other CMASS measurements, even 
though our theoretical models are fairly different. This consistency suggests that measurement 
of fag from Redshift space distortions at multiple redshifts will be a sensitive probe of the 
theory of gravity that is largely model independent, allowing us to place model-independent 
constraints on alternative models of gravity. 

Key words: gravitation; galaxies: haloes; galaxies: statistics; cosmological parameters; large- 
scale structure of Universe 


1 INTRODUCTION 

The evolution of our Universe ap pears to be wel l described by the 
theory of general relativity(GR) jPeeblesI I198QI : IPavis & Peebles! 
Il983h . The predictions appear to be consistent with all the obser- 


vations except the mysterious accelerated exi 

mansion of the Uni- 
b, and the dark mat- 

verse (Riess et al. 

19981; Perlmutter et al.lll99^ 

ter JZwickvl 1937 

;lKahn&Wolfiei]ll959l; Rubin & Pordlll970h . 

The accelerated expansion of the Universe proposed in ACDM- 
GR is in good agreement with Baryon Acoustic Oscillation (BAO) 

(Eisenstein et al.1 20051; Cole et al.l 

20051; Hutsill200fil: iKazin et all 

I 2 OIOI: IPercival et al.l I 2 OIOI; Reid et al. I 2 OIOI 

), Hubble constant 
Background (CMB) 

jRiess et al.l 201 lb and Cosmic Microwave 1 


teennett et al .120131 : IPlanck Collaboration et al.l2014ah . Within the 


current paradigm, the primordial fluctuations in the early Universe 
were amplified into structures we observe today via gravitational 


interactions. These gravitational interactions are the sum of the 
motions of matter and the expansion of the Universe. Therefore, 
one would only need to specify the initial conditions, the spatial 
geometry and the contents of the Universe to use Einstein’s the¬ 
ory to predict the large-scale growth rate of the matter density in 
the Universe. We can co mpare such predictions to the observa¬ 
tions in redshift su r veys dSchlegel et'Sl l2009l: iBlake et a , |l201ll: 
iBeutler et alJ 1201 2 I : IPercival et al. 1200-1 Ide la Torre et al. 2013h . 
The velocity field from maps of galaxies in such surveys can be 
measured because the galaxy redshifts, from which distances are 
inferred, include components from both the Hubble flow and pe¬ 
culiar velocities from the comoving motions of galaxies . Such sur¬ 
veys thus reveal an anisotropic distribution of o bjects (ICole et al.l 
ll995l : IP^cock et al1l200ll : [^mushia et~^l201lh : the anisotropy in 
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the clustering encodes information about the formation of structure 
and provides a sharp test of the theory. 

In galaxy redshift surveys, the distortion produced in the two- 
point correlation function due to the peculiar velocity compo¬ 
nent in the g alaxy r e dshift is known as “Redshift Space Disto- 
rion (RSD)”. iKaised (Il987h first developed a formalism that de¬ 
scribes redshift space power s pectrum by modi fying the linear the¬ 


ory of large scale structure. iHamiltonf (Il992h extended this ap¬ 
proach to the tw o-point correlation f unction in real space. The 
seminal work of IScoccim^ 3 (l2004 describes a more general 
dispersion model, which imp roves the Ka i ser lin ear model by 
including higher order terms. IScoccima^ (l2004l) also invoked 
the conc ept of general st r eamin g model which was first intro¬ 
duced by IPavis & PeeblesI (Il983h and further developed by iFishej 
(Il995h . A combination of th e Lagrangian Perturbation Theory 
model by IScoccim^^ (l2004h and the Gaussian Streaming model 
was used to measure the linear growth rate (/erg) o f the Universe 
in Ba ryon Oscillation Spectroscopic Survey fBOSS: [Schlegel et~^ 
l2009h. which is part of Sloan Digital Sky Survey I II (SPSS- 
III: lEisenstein et al.ll201 ih. DR9 dReid & WhitelEonI: iReid et alJ 
l2012h and DRl 1 dSamushia et alJlioTsl) data releases. Other meth- 
ods and models of the correlation functions or power speetruin 
are also used to der i ve fag with BOSS d ata (IChuang et al.ll2013] : 

2013 ). Other galaxy redshift 
' ^,6-d egree Field 
T 2012 I)). 2-degre e 


Ib cutler etaP l20l4 ISanchez et al 


ett et al.l j20 


Bcutler et al. 


surveys such as SPSS DR7 iHow 

Galaxy Redshift Survey fhdFGRS _ _ - 

Field Gal axy Redshift Surv ey r2dFGRs lPercival et al.l (l2004l) ) 
WiggleZ telake et al.l 12011 ), and VIMOS Pub lic Extragalactic 


Redshift Survey fVIPERS de la Torre et al.l (l2013h ) have also mea¬ 
sured redshift space distortion at different redshifts. 

In this paper we employ a model called Convolution La¬ 
grangian Perturbation Theory with Gaussian Streaming Redshift 
Space Distortions (hereafter CLPT-GSRSD) to analyze BOSS 
DRll (I Anderson et al.l 12014^ . In order to test the model and 
the analysis method presented in this paper to a high accuracy, 
we have used a rela tively large number of mock galaxy catalogs 
jManera et al .1120131) with clustering properties similar to those of 
the higher redshift BOSS galaxies. We also provide systematic er¬ 
rors based on the results of analyzing these mock galaxy catalogs. 


2 THEORY 

In this section we describe the model with which we fit both the 
mock galaxy catalogs and BOSS data. 

Throughout the paper, we adopt the standard “plane-parallel” 
or “distant-observer” approximation, in which the line-of-sight di¬ 
rection to each object is taken to be the fixed direction i. This ap¬ 
proach has been shown to be a good approximatio n at the level of 
curren t observational error bars (e-g-, Figu re 10 of ISamushia et al.l 
(l2012h or Figure 8 of lYoo & SeliakrfcOlSh ). 


2.1 CLPT 

Convolution Lagrangian Perturbation Theory (CLPT) is a non- 

i )erturbative resumm ation of Lagrangian perturbation theory 
Carlson et~aDl2013h . With CLPT, Carlson identified a few terms 
that asymptote to constants in the large-scale limit and hence need 
not be expanded with approximations. The authors showed that 
CLPT performs better than all the other methods when compared 
to N-body simulations of dark matter halos. The monopole of cor¬ 
relation function matches N-body up to a very small scale; the 


quadrupole has less than a few percent error f or sca le above 20 
h“^Mpc (see Figure 1 and 2 of Icarlson et al.l (l2013h ) . Unfortu¬ 
nately, the CLPT doesn’t perform well in the quasi -linear regime 
for the quadrupole of biased tracer (see Figure 5 of ICarlson et al.l 
(l2013h ). To overcome this problem with the quadrupole for a biased 
tracer we use CLPT in combination with the Gaussian Streaming 
Model (GSM) as it has been demon strated to model th e galaxy cor¬ 
relation functions more accurately dWang et al.ll2014h . 


2.2 The Gaussian streaming model 


The Gaussia n Str eaming Model (GSM) developed by 
iReid & White! (l201ll) . fits the monopole and quadrupole of 
the correlation functions of mock galaxies with a large-scale bias 
b ^ 2io the percent level on scales above 25 /i“^Mpc. This model 
has be en used to interpret the clustering of galaxie s measured in 
BOSS dReid et al.ll2014ISamushia et al.l2013ri2014h . 

The GSM is inspired by the Eulerian streaming models. It en¬ 
forces pair conservation, assuming that the functional form of the 
halo velocity distribution is Gaussian, centered at 11 V 12 , where //U 12 
is mean line of sight velocity between a pair of tracers as a func¬ 
tion of their real space separation. Specifically we assume that the 
redshift-space halo correlation function is 


l+^^(s±,S||) = 


-I 


dy 


V^ai2 


[l+^(r)]exp< 


[g|l - y- yvi 2 ]‘^ \ 

+ ^fog) J 
( 1 ) 


where ^(r), V 12 and a 12 are provi d ed from an analyt ic theory. In 
the model of iReid & White! d201lh : iReid et all d2012h intergrated 
Lagrangian pe rturbation theory with scale-dependent but local La¬ 
grangian bias dMatsubarall2008l) was used for the real-space corre¬ 
lation function (^(r)) of halos, but the halo infall velocity (vi 2 {r)) 
and dispersion(cri 2 (r)) were computed in standard perturbation 
theory with scale-independen t bias. In order to move from halos 
to galaxies, iReid et all d2012h showed that it suffices to introduce 
a single additional parameter, a fog, akin to the a in Eq. 0. 
This quantity is taken to be an isotropic, scale-independent dis¬ 
persion that is added in quadrature to a 12 that modifies the scale- 
dependence of the quadrupole moment on small scales. There 
are more comprehensive simulation based models for describing 
the velocity distributi on of galaxies aro und groups and clu sters 
dZu & Weinberg||2013h and at large scales dBianchi et al.l2015h . but 
it is relatively difficult to embed them into the halo model to explain 
the kinematics of the galaxies. 


2.3 CLPT-GSRSD 

The quadrupole prediction on quasi-linear scales from CLPT for a 
biased tracer is not as accurate as its predictions for N-body sim- 
ulations. CLPT w as further improved by the model proposed by 
IWang et ^ d2014h . which combines the velocity statistics and cor¬ 
relation function from CLPT with GSM in order to produce a more 
accurate monopole and quadrupole for biased tracer. This model 
has less than 5% error in quadrup ole for pair separ ation greater 
than 20 h“^Mpc for a biased tracer dWang et al.ll2014l). This mode l 
is similar to the Zeldvoich Streaming Mo del (ZSM) ( Whitel[20l4) 
and Lagrangian Streaming Model (LSM) dWhite et al.l2015h . None 
of these models has been yet used to interpret the clustering of the 
BOSS gal axy sample. We w ill be using the analytical model de¬ 
scribed bv lWang et al.l d2014l) to extract constraints on cosmologi¬ 
cal parameters; we denote this model CLPT-GSRSD. 
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3 SDSS III- BOSS DATA 


We use data i ncluded in data rele ases 10 fPRlO lAhn et alJ (l20l40 
and 11 tPRlllAlam et alJ (l2015h ) of the Sloan Pi gital Sky Survey 
(SPS S; lYork et alJ (l200Qh\ Together, SPSS I, II (lAbazaiian et alJ 


l2009h and III jEisenstein et al. 201l[) used a drift-scanning mo¬ 
saic CCP camera (iGunn et alJ 1998h to image over one-third of 


the sky (14555 square degrees) in five photometric b andpasses 
jFukugita et al.ll 19961 : ISmith et ^l2002l : |Poi et al.ll201Qh to a limit¬ 
ing m agnitude of r < 2 2.5 using the dedicated 2.5-m Sloan Tele¬ 
scope ( Gunn et aDl2006h located at the Apache Point Observatory 
in New Mexico. The imaging data were processe d through a serie s 
of pipelines that perfor m astrometric calibr ation jPier et al.ll2003h . 
photom etric reduction (iLupton et al.lll999h . and photometric cali¬ 
bration jPadmanabhan et al .1120081) . All of th e imaging was repr o- 
cessed as part of SPSS Pata Release 8 fPR8: lAihara et al.l (l201lh ). 
BOSS (jPawson et al.ll201!^ is designed to obtain spectra and red- 
shifts for 1.35 million galaxies over a footprint covering 10,000 
square degrees. These galaxies are selected from the SPSS PR8 
imaging and are being observed together with 160,000 quasars and 
approximately 100,000 ancillary targets. The targets are assigned 
to tiles using a tiling algorithm that is adaptive to the density of tar¬ 
gets on the sky (iBlanton et al.l2003]). Spectra are obtai ned using the 
double-armed BOSS spectrographs (ISmee et al.l2013h . Each obser¬ 
vation is performed in a series of 900-second exposures, integrat¬ 
ing until a minimum signal-to-noise ratio is achieved for the faint 
galaxy targets. This ensures a homogeneous data set with a high 
redshift completeness of more than 97% over the full survey foot¬ 
print. Redsh ifts are extracted fr om the spectra using the methods 
described i niBolton et al. (|2012h. A summary of the survey design 
appears in Eisenstein et al.l~( 201 ih . and a full description is pro¬ 
vided in IPawson et al. (12013 ). 


We use the CM ASS sample of ga laxies dBolton et al.ll2012h 
from data release 11 (lAlam et al .112015h . The CMASS sample has 
690,826 massive galaxies covering 8498 square degrees in the red- 
shift range 0.43 < z < 0.70, which correspond to an effective 
volume of 6 Gpc^. 


4 INVESTIGATING THE SYSTEMATICS BUDGET 

In this section we will examine possible sources for the systematics. 
We will first describe the result of fitting PTHalo mocks and then 
the systematic introduced by various approximations used in the 
analysis. We will conclude this section by describing the possible 
observational systematics. 


4.1 Mock Galaxy catalogs 

To validate CLPT-GSRSP we have use d PTHalo niock ca talogs. 
Such validations have been performed bv lWang et al.l (l2014l) in the 
paper in which they were introduced. However, our goal is to test 
the model for mock galaxies with properties similar to those of the 
BOSS galaxies at z ^ 0.5. 

The perturbation theory mock (PTHalo mock) is generated by 
populating a matter field using second order Lagrangian pertur¬ 
bation theory and calibrating the masses of da rk matter halos by 
compying it to detailed n umerical simulations (IScoccimarro et al.l 
l200lHManeraetalJl2013h. We are u s ing th e PTHalo mocks gen¬ 
erated and validated bv iManera et al.l (l2013h for PRl 1 footprint of 
SDSS-III survey. There are 600 PTHalo mocks available, which we 
employed to validate the model and decide upon some constraints 


on the model parameters. The simulations used to produce PTHalo 
mocks cover the same volume as that of the CMASS sample. These 
mocks are designed to have a bias similar to the bias of CMASS 
sample. It is important to note that both, mocks and our model, are 
based on a similar kind of perturbation theory. 

4.2 Fitting PTHalo mocks 

We have examined the result of fitting 600 PTHalo mocks. We 
fixed the cosmology while fitting these mocks because we pre¬ 
cisely know the cosmology of simulation. The mocks are con¬ 
structed in such a way that they mimic the CMASS sample . We 
have computed the monopole and quadrupole of the correlation 
fun ction for each of the 6 00 mocks using Landy-Szalay estima¬ 
tor dUandv & Szalavl[r993h with bins of 8 h“^Mpc in pair separa¬ 
tion. The monopole and quadrupole are fit with the CLPT-GSRSD 
model for fixed fiducial cosmology and three freely-floating RSD 
parameters, {F', F", /}, where the parameters represent the first 
and the second order Lagrangian bias and the derivative of loga¬ 
rithm of the growth factor, respectively. Figure[T]shows the result of 
fitting 600 PTHalo mocks using CLPT-GSRSD model as the func¬ 
tion of minimum scale (smin) used for the fit. The top panel demon¬ 
strates that we recover the expected value of / = d{\n D) /d{\n a) 
within a few percent for the entire range of Smin shown in the fig¬ 
ure. The bottom panel reveals that the /do f increases as we in¬ 
clude smaller scales in the fit. The solid blue lines in these plots 
are the mode of the results from 600 mocks and the shaded region 
corresponds to la limit. We have also performed this analysis with 
bins of 2 h~^Mpc in mon opole and quadrupole and found similar 
results. IWang et ^ (l2014l) has shown that the prediction of corre¬ 
lation function by CLPT-GSRSD agrees at few percent level down 
to r = 15 h“^Mpc with N-body simulation. We would like to keep 
Smin as small as possible because the correlation function has a 
high signal to noise at smaller scales. But, our mocks are based on a 
perturbation theory which means the shape of quadrupole might not 
be very accurate at small scale due to the non-linear velocity struc¬ 
ture. It is clear from Figure [T] that below 36 h“^Mpc the growth 
rate estimated starts to get slightly biased and /dof crosses 1. 
So, this is similar to bias-variance trade off, where Smin = 36 h“^ 
Mpc shows no bias and Smin = 22 h“^Mpc is with smaller vari¬ 
ance. We have decided to use Smin = 30 h“^ Mpc as the optimum 
minimum scale with acceptable bias and low variance in our anal¬ 
ysis. We have also repeated the analysis with Smin = 22 h“^Mpc 
and found consistent result with marginal improvement in preci¬ 
sion. 

Figure |2]demonstrates that the observed scatter in the parame¬ 
ter produced by fitting the mocks is completely consistent with the 
model degeneracy. The black points are the best fit obtained from 
the mock (some of the mocks are outside of the region shown in 
the plot). The red dashed contours in the figure are an estimate of 
the model degeneracy. We have evaluated the theoretical correla¬ 
tion function on the 100 x 100 grid in the region covered in Figure 
|2] We have used correlation function evaluated for F' = 1.0 and 
/ = 0.76 as the reference model to evaluate the ^^ch point 
of the grid using the covariance matrix evaluated from the mocks. 
The x2 surface constructed in this manner should reveal the de¬ 
generacy of the model. The red dashed line in Figured represents 
the contour of x^ while using fitting scale 30 h“^Mpc < s < 126 
h“^Mpc. We designate this as model degeneracy and show that 
almost all the scatter in the mocks can be explained by the model 
degeneracy. We have found 17 mocks with Ff < 0.8 (outliers). The 
mean /do f for the outlier mocks is 7.7, which is relatively high. 
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Figure 1. The result of fitting 600 PTHalo mocks. The two panels show the 
mode and la (computed from percentile ) spread of {/, x'^ldof} as the 
function of Smin minimum scale used for fitting the correlation function. 
In the top panel for /, we recover the expected value of the parameter within 
a few percent. The vertical dashed line indicates the minimum scale used 
in this paper Smin = 30 h“^Mpc. The lower panel demonstrates that the 
model and data do not agree as we move to smaller(more non-linear) scale 



Table 1. The shift in the parameter {fas, Q'||, for each weight ap¬ 
plication from the default of applying all weights in units of error on the 
parameter. The stellar weight is most important for the full shape analysis. 


Weights 

A/(T8/(T/<jg 

Aq;|| /crQ:|| 

AQ;^/crc_L 

'^star * 1) 

-0.03 

0.02 

0.03 

Wsee * {Wcp + Wzf - 1) 

0.54 

0.82 

0.12 

Wcp + Wzf - 1 

0.51 

0.82 

0.13 

'^zf 

0.79 

0.71 

-0.20 

Wcp 

0.44 

0.94 

0.36 

None (uniform weighting) 

0.72 

0.82 

0.06 


the weights to account for redshift failure (Wzf) and close pairs 
(Wcp). Wcp accounts for the fact that the spectrum of only one of 
the galaxies is observed for galaxi es separated by les s than 62” due 
to the finite size of fiber fittings. IRoss et aP (l2012h demonstrates 
that these weights don’t change the cosmological result for B ary on 
Acoustic Oscillation (BAO) but do remove the systemati c effects in 
the me asurement of and improve the goodness of fit. IRoss et ahl 
(l2012h have also shown that for scales above 150 h“^Mpc has 
a systematic error larger than the statistical error. This issue is not 
a concern for our study because our maximum scale of interest is 
126 h“^Mpc. 

We have also examined the effect of usin g different w e ightin g 
on the measured RSD parameters (similar to lOsumi et al] (l2015h ) 
using observed correlation function. We have calculated the corre¬ 
lation function with different weighting, these correlation functions 
are fit for RSD and BAO parameters [fag, bias, apoG, tty, <^±}. 
Table [T] lists the observed shift compared to the default case of us¬ 
ing all weights. The stellar weight has the important effects in our 
measurements. If we do not use any weights the shifts are of the 
order of la for fag and ay but is small for a±. The shifts in mea¬ 
surement due to these systematic weights will be higher as we will 
improve precision in future surveys. 


5 ANALYSIS 


Figure 2. The scatter shown by black points is the result of fitting 600 
PTHalo mocks. Seven out of 600 mocks showed a best fit that is outside the 
region shown in the plot. The red-dashed contour lines represent the model 
degeneracy. This plot demonstrates that most of scattering in the mock can 
be explained by the model degeneracy. We have found 17 outlier mocks 
with F/ < 0.8 and relatively high x^- 

The mean of quadrupole moment of these mocks shows stronger 
variation around r = 100 h“^Mpc. 

4.3 Observational effects 

There are several possible sources of observational systematics in 
the measurement of correlation function. These systematics have 
be en studied i n grea t detail for the SPSS DR 9 and DR 11 samples 
in IRoss et^ (l2012h : lAnderson et al.l (l2014bl l3). They have intro¬ 
duced two different systematic weights to reduce the cross corre¬ 
lation of with the star density and seeing, while demonstrating 
that potential systematics, such as sky brightness, do not affect the 
measurements. We have adopted the weights Wstar and Wsee in our 
measurement to remove these systematic effects. We have also used 


In this section we describe the steps of our analysis starting from 
galaxy position to the parameter constraints. We first briefly men¬ 
tion the different steps then elaborate on the procedure with details 
in the following subsections . 

We first measure the galaxy correlation function using SDSS- 
DRll (CMASS) galaxy sample. We have used 600 PTHalo mocks 
to generate an estimate of the covariance matrix. We start our opti¬ 
mization problem with nine-dimensional parameter space consist¬ 
ing of four cosmological parameters and three RSD parameters, 
and run Markov Chain Monte Ca rlo(MCMC) to explor e this pa¬ 
rameter space using COSMOMC (iLewis & Bridlel[2003) . In every 
MCMC step, we first evaluate the Planck likelihood, and then eval¬ 
uate the linear power spectrum for the current sampled(grid) point 
if it has not been evaluated previously. This linear power spectrum 
is fed to the CLPT theory to calculate the correlation function(^(r)) 
and velocity statistics. The CLPT ^(r) and velocity statistics are 
used by the Gaussian streaming model GSRSD with the three RSD 
parameters to evaluate the redshift space two-dimensional correla¬ 
tion function. This two-dimensional correlation function is rescaled 
according to the difference in the fiducial cosmology and the cur¬ 
rent MCMC cosmology to determine the final model correlation 
function. The final correlation function is used as the theoretical 
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model (CLPT-GSRSD) to calculate the Xrsd with data correla¬ 
tion function and mock covariance matrix. The Xrsd is converted 
to a likelihood, which is multiplied to the Planck likelihood in order 
to calculate the total likelihood, which is maximized using COS- 
MOMC. 

The basic approach is to estimate the cosmological parame¬ 
ters by fitting the monopole and quadrupole of correlation func¬ 
tion from SDSS CMASS DRll using the correlation function of 
CLPT-GSRSD in combination with Planck likelihood computed 
from CMB power spectrum. 


5.1 Measuring the correlation function and covariance 
matrix 

We first assumed the fiducial cosmology as flat ACDM-GR cos¬ 
mological model with Qm = 0.274, Hp = 0.7 ,p bh^ = 0.0224, 

Us = 0.95 and ag = 0.8 (I Anderson et al.ll^l4ah in order to con¬ 
vert observed celestial coordinates (a, S) and redshift to the posi¬ 
tion of the galaxy in three-dimensional space. These galaxy posi¬ 
tions are used to estimate the two point statistic (correlation func¬ 
tion) of the galax y using the minimum variance estimator Landy- 
Szalay estimator (ILandv & Szalavl[l9^ . 

^ ^ DD{AT)-2DR{Ar)^RR{AT) 

^ RR{Ar) 

where DD, DR and RR represent weighted pair count of galaxy, 
cross pair count of galaxy-random and pair counts of randoms, re¬ 
spectively. We used the weighting wtot = WfkryW starWseeiwcv + 

Wzf — 1) as described in lAnderson et al.l(l2014ah to correct for ob¬ 
servational systematic errors and minimize variance. The correla¬ 
tion function is function of r, which is the distance between a pair 
of galaxies and ji — cos{0), where 0 is the angle between separa¬ 
tion vector and line of sight. We compress the information by pro¬ 
jecting the correlation function to the Legendre polynomial Li (/x) 
of order I as follows. 

iiir) = ‘^^2^ [ dni{r,n)Li{n) fa?^^^'^Anki{r,^i)Li{nk) 

(3) 

IWang et aP (I 2 OI 4 I) reported that the CLPT-GSRSD model is 
a good fit to N-body simulation for / = 0 and I = 2 but not for 
higher /, therefore, we will limit our analysis up to / = 2. The 
signal measured in the correlation function should evolve with red- 
shift due to the evolution of osiz), but a good approximation is to 
use an effective re dshift for DRll sampl e z = 0.57 as shown in 
equation[10-12] of ISamushia et al.l (l2013h . 

To estimate the uncertainty in our measurement of the cor¬ 
relation function, we constructed the covariance matrix using 
600 PTHalo mocks. The inverse of the sample covariance esti¬ 
mated from finite number of mocks is a biased estimator. We 
will adopt the corre c tion f o r covariance m atrix as described in 
IVargas Magana et aP (l2013h : IFercivall (l2013h . The sample covari¬ 
ance matrix is calculated as follows. 


A m - 0 + El=30l(Ci - m - 0 

2 X 299 

(4) 

where Cij represents the covariance between bin i and j and ^ is 
the mean of the mocks and the sum is over different mocks. The 
corrected inverse covariance matrix is 


c~j = (1 - 0.62(2 X nr„„ + l)/{N^ock - (5) 

where Nmock is the total number of mocks which is 600 for us, and 
Tiru is the number of bins of the correlation function used in the 

' bin 

analysis. 


5.2 Parameter space 

Our model parameters can be subdivided into two subsets. The first 
set is of cosmological parameter ns,As Hq), where 

Ob, Oc and Ho are baryon density, dark matter density and bub¬ 
ble constant, respectively, with h being Ho/100. The quantity As 
is the scalar amplitude of primordial power spectrum. This choice 
of parameters requires us to assume an evolution model in order to 
evaluate the theoretical model at the galaxy redshift.We want our 
measurement to be independent of such an assumption, which is 
necessary to be able to use these results to test various models of 
gravity. This goal can be accomplished by fixing Ho to the best fit 
and allowing two extra parameters, H{z) and Da{z), at the effec¬ 
tive redshift. This can be modeled using Alcock-Paczynski param¬ 
eters The second set of parameters are the redshift space 

distortion (RSD) parameters {F',F" J.gfog), where F\F" and 
/ represent the first and the second order Lagrangian bias, and the 
logarithmic derivative of the growth factor, respectively. The pa¬ 
rameter a fog is an additional is otropic velocity disp ersion to ac¬ 
count for the finger-of-god effect (iReid & Whitel201 u) . The second 
order Lagrangian bias is not well constrained because it is more im¬ 
portant for small scale quadrupole. We have considered two cases. 
In the first one we have marginalized over F" with hard prior cov¬ 
ering [—5, 5] in the second case we have sampled overdensity r 
which determ ines both F' and F" using peak background split 
(IWhitell20l4 . 


F' = 
F" = 


Tc 


2 1 , 
av — 1 + 


1 + {aFy 


2 4 o 2 , 2p(2aF + 2p - 1) 
a R — 3aiy H-^^ 


1 + 


( 6 ) 

(7) 


where a = 0.707, p = 0 .3 gives the Sheth-Tormen mass function 
(ISheth & Tormenlfl^^ . and 5c = 1.686 is the critical density for 
collapse. Independent of whether we use peak background split or 
marginalize over F", we obtain consistent results . 


5.3 Calculating Theoretical model 

We first calculate the linear matter power spectrum usi ng Code for 
Aniso tropies in the Microwave Background (CAMB) (iLewis et al.l 
l200Qh . This approach requires the knowledge of cosmological pa¬ 
rameters of the model we are evaluating. The linear power spectrum 
is then sampled according to the sampling scheme for k described 
in Appendix IaI The sampled power spectra is used to calculate the 
real space correlation function (^(r)), pairwise velocity statistics 
(vi 2 {r)) and the dispersion of pair wise velocity statistics ( a i 2 (r) 
) as described in IWang et al.l (l2014l) and ICarlson et al.l (l2013h . 
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'^ + £.{r) = j d^qMo{r, q) (8) 

vi2(r) = [1 + CW]~^ J d^qMi,n{r,q) (9) 

0-12,nmW = [l+CW]~^y d^qM2,nm{r,q) (10) 

= '^(^12,nrarnTm (H) 

nm 

<^l{r) = ^((T?2,nm5nm “ C^y)/2 (12) 

nm 


where Mo{r,q) and M2^nm{r,q) are integrals of the 

CLPT perturbation theory that depend on the linear matter power 
spectra. The fn ,fm are unit vec t ors alo ng the galaxies position vec¬ 
tors. Please refer to I Wang et al.l (l2014l) for details of these integrals 
and derivation of these equations. 

The real space correlation function and velocity statistics cal¬ 
culated from CLPT together with the RSD parameters are suppl ied 
to the Gaussian Streaming Model (GSM) dReid & Whitell201 ih in 
order to evaluate the redshift space correlation function as follows. 


i+r°‘''^'(s±,sii) = 

f dy[l + ^{r)] .f [S||-ry-Mt’i 2 (r)]" 

J [27rcr22(r,^)]i/2 P\ 2cr?2(r,^) 

a’i2(?',M) = W + (1 - + <^FOG (14) 

where s± is the transverse separation in both real and redshift 
space. The quantities s\\ and ry are the LOS (line of sight) sepa¬ 
ration in redshift and real space, respectively. 

We want our measurement to be independent of any particular 
model of gravity. This can be achieved by avoiding the use of any 
model for the evolution of structure formation under the assump¬ 
tion of small deviation from widely accepted ACDM-GR. We will 
model this deviation by using two parameters (ay, a±) which are 
defined as follows. 



dimensional parameter space. We have four cosmological param¬ 
eters {Qbh^^ Qch^, Us, As}. These parameters have flat prior cen¬ 
tered at best fit value of Planck with width :£5apianck (see Table 
|2]). The other set of parameters are five redshift space distortion pa¬ 
rameters {ly, /, a fog, Q^ii , where u is overdensity, which de¬ 

termines the first and the second order bias using peak background 
split relation. When we marginalize over second order Lagrangian 
bias F" in place of using the peak background split, then our sec¬ 
ond set of parameters is replaced by {F' f,aFOG,G\\,a±}. 

We calculate the model monopole and quadrupole of the 
galaxy correlation function as described in the previous section for 
each point in the parameter space visited by the MCMC sampler. 
We have optimized our model evaluations with some assumptions, 
discussed in the Appendix |3 in order to make time per MCMC 
iteration smaller. The likelihood constraint from RSD is calculated 
using as follows. 


XrSD — {^model ^data) C {^model ^data ) (17) 

Crsd = exp(-x^/2) (18) 

where ^m.odei = , Uata = and C~^ 

is the inverse of corrected covariance matrix calculated from 600 
PTHalo mocks. The cosmological parameters are well constrained 
from Planck satellite CMB data, which are mostly independent of 
the gravity and growth of structure parameters. Therefore, we will 
multiply our RSD likelihood by the Planck likelihood to obtain the 
joint constraint on our parameters. 


^total — Fplanck^RS D (19) 

The Planck likelihood is calculated using the constraint 
reported from planck t e mperatu re anisotropy data alone 
(jPlanck Collaboration et all 12014bh . The Planck parameter 
covari ance is obtained from the correlati on matrix given in Figure 
21 of jPlanck Collaboration et al.ll2014ch . We employ a multivari¬ 
ate Gaussian approximation to the full Planck likelihood. This 
approximation is close to the actual likelihood in the parameter 
space in which we are working. 


^fiducial /i c\ 

= - ^,«x= Xdull (15) 


Hf 




D 


We will rescale the model redshift space correlation function 
to account for this extra distortion as follow. 


= 

l-s —So I < As /2 


-Number of bins used in sum 

TT 


(16) 


where {Hfiducial , H} and ^ are the bubble expan¬ 

sion rate and angular diameter distance for the fiducial and model 
cosmology, respectively, So = As = 5 

h“^Mpc. The ab ove rescaling is simply the application of Alcock- 
Paczynski effect (lAlcock & PaczvnsS|[l979h . 


5.4 MCMC and Likelihood estimation 

We use COSMOM C to perform Markov Chain Monte-Carlo 
likelihood analysis (iLewis & Bridlel l2002h and explore nine- 


Vtbh^ = 0.02207, Vtch^ = 0.1196, n, = 0.9616, A, = 3.098 


Cplanck — 


/ 1.089X10“'^ 
-4.501X10“'^ 
1.365X10“® 

\ 3.564x10“® 


-4.501X10“'^ 1.365X10“® 3.564x10“® \ 

9.610X10“® -2.215X10“® 1.562x10“® 

-2.215X10“® 8.836X10“® 2.030x10“® 

1.562x10“® 2.030X10“® 5.184x10“®/ 


6 RESULTS 

Figure [3 presents the monopole and quadrupole of the galaxy cor¬ 
relation function. The black data points are measurements from 
BOSS CMASS DRll galaxy sample. The error bars on the mea¬ 
surements are the diagonal elements of mock covariance ma¬ 
trix. The blue and red line represents the best fit monopole and 
quadrupole, using the fitting range 30 h“^Mpc ^ s ^ 126 
h“^Mpc with 8 h“^Mpc sampling. The best fit /do f = 24/22 
is achieved after marginalizing over essentially all the relevant pa¬ 
rameters as listed in Table |2] The Figures |4] and [5] show the one di¬ 
mensional and two dimensional marginalized likelihood for some 
of the parameters. The final results of this analysis are given in 
Table |2] We have measured fasiz = 0.57) = 0.462 zb 0.041 
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Table 2. The list of parameters used in our analysis. For each parameter we provide their symbol, prior range, central value and Icr error. We report the results 
of using peak background split to relate the first and second order bias and also the result of using sampling of the two bias independently. The results in both 
case are consistent. We also list the result when Planck prior is replaced by WMAP prior, which predicts 2% shift in fag. 


Parameter prior range Peak background split First and Second order bias 

with WMAP prior with Planck prior with Planck prior 


Sampling Parameters 


Ubh'^ ... 

[0.02042,0.02372 ] 

0.02267 ± 0.00036 

0.02206 ± 0.00026 

0.02206 ± 0.00026 

... 

[0.1041,0.1351] 

0.1141 ±0.0021 

0.11956 ±0.00086 

0.11956 ±0.00086 

ris ... 

[0.914,1.008] 

0.9741 ± 0.0085 

0.9614 ± 0.0058 

0.9613 ± 0.0058 

ln(10“As)... 

[2.67,3.535] 

3.178 ±0.029 

3.093 ± 0.066 

3.103 ±0.070 

ay 

[0.8,1.2] 

1.003 ± 0.039 

1.051 ±0.043 

1.058 ±0.047 

Ox 

[0.8 ,1.2] 

0.997 ±0.018 

1.03 ±0.016 

1.032 ±0.016 

/ = dlnD / dlna 

[0.3 , 1.2] 

0.739 ± 0.067 

0.747 ±0.072 

0.729 ±0.073 

C^FOG 

[0,10] 

2.26 ± 1.46 

1.91 ± 1.28 

2.70 ± 1.69 

vrsd 

[1.5,2.0] 

1.83 ±0.038 

1.80 ±0.05 


F' 

[0.5 , 1.5] 



0.93 ±0.07 

F" 

[0.5 , 1.5] 



1.0 ± 1.78 

Derived Parameters 

/cts 


0.454 ± 0.041 

0.462 ± 0.041 

0.453 ± 0.041 

bus 


1.21 ±0.030 

1.194 ±0.032 

1.20 ±0.032 

= 0.57) 


1356.0 ±24.0 

1400.9 ± 22.7 

1403 ±21.9 

H(z = 0.57) 


93.4 ±3.6 

89.2 ±3.6 

88.5 ±3.9 

Fap 


0.663 ± 0.033 

0.654 ± 0.033 

0.651 ±0.034 

Dv{z = 0.57) 


2024.5 ± 27.2 

2101.4 ±25.6 

2108 ±29.4 




29 0.40 0.51 0.62 




Figure 3. This hgure shows the projection of two-dimensional correlation 
function in legendre basis. The black data points are measured correlation 
function for DRl 1 CM ASS sample and the error bars are the diagonal terms 
of covariance matrix calculated using 600 PTHALO mocks. The blue and 
red line represents the best ht monopole and quadrupole with ^ jdo f = 
24/22. 


,6(78 = 1.19 zb 0.03, Da{z = 0.57) = 1401 zb 23 Mpc and 
H{^z — 0.57) = 89.2 zb 3.6 km Mpc“^. The galaxy correla¬ 
tion function doesn’t improve the constraints on the baryon density 
scalar spectral index(ns) and amplitude of primordial cur¬ 
vature perturbation (As) at ko = 0.05 Mpc“^h over the already 
tight constraints from Planck, however we do improve the mea¬ 
surement of cold dark matter density (Qch^ = 0.1196 zb 0.0009) 
,which is a 0.7% measurement, this is an improvement in measure¬ 
ment of cold dark matter density compared to Planck (Qch^ = 
0.1196 zb 0.0031) by a factor of 3.6 . We have also studied the im¬ 
pact of second order bias on growth rate measurement. We found 
consistent constraint while allowing both of first and second or¬ 
der Lagrangian bias to be free (shown by blue dashed line in Fig- 


1295.0 1364.9 1434.8 1504.7 

Da (^=0.57) 

— Planck+peak 
—- Planck+F1,F2 

Figure 4. The four panels show the one-dimensional marginalized likeli¬ 
hood for the parameters {fag^bag, Da, H} at effective redshift 0.57 . 
The red solid line shows the result with peak background split, and the 
blue dashed line is the result when we fit for both first and second order 
Lagrangian bias. The grey shaded region shows Icr constraint from Planck 
with ACDM-GR. We detect 1.9% shift in fag and less than 1% for other 
parameters between the analysis with and without peak background split. 


ure 0. It is interesting to note that we have found covariance be¬ 
tween second order Lagrangian bias (F") and growth rate. This 
suggest that a better modeling of higher-order bias including local 
and non-local contrib utions (for a recent work along this line, see 
e.g. lSaito et al.l (l2014l) ) will become more important as the cluster¬ 
ing statistics become more precise with future surveys.We have also 
repeat ed our analysis by r eplacing the Planck prior with WMAP 
prior teennett et al.l 120131) and found 2% shift in fag, which is 
much smaller than the estimated error. 


1.0 ..34 



73.37 84.11 94.85 105.59 

H{z =0.57) 
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i.OOr^ , , , 

0.25 0.38 0.52 0.65 



13lb'''' 1367'" 1423 ■ " i'480 


Da (^ = 0.57) 


Figure 5. Two-dimensional 68% (la) and 95% (2a) confidence limits ob¬ 
tained on fag—bag and Da—H at effective redshift of 0.57 recovered 
from Planck CMB and CM ASS (^ 0 , 2 ) datasets with peak background split 
assumption. 


0.55 



u 


0.30^-^ 

Figure 6. Comparison of fag with other analysis on the same DRll 
CMASS sample. The blue point present the result from our analysis. 
Our measurement is consistent with other clustering analysis and Planck 
ACDM-GR prediction. 


Our measurement of fag is consistent with all the other mea- 
surements repo r ted fro m the same data set, as shown in Figure [b] 
ISamushia"^ "S] (l2013h has reported fag{z = 0.57) = 0.44 zb 
0.044, which is similar to our analysis. One major difference is 
in the theoretical model used in the two studies. We have used 
CLPT-GSRSD t o evalu ate our model correlation function, whereas 
ISamushia et al.l (l2013h use Lagrangian pertu rbation theory (LPT) 
as the model to predict correlation function. Ib cutler et al.l (l2014h 
used the monopole and quadrupole of power s pectrum and reported 
fag{z = 0.57) = 0.419 zb 0.044, whereas Ichuang et ^ J2013h 
performs the analysis in configuration space with a different fit¬ 
ting m odel and obtained fag = 0.391 zb 0.044. ISanchez et al.l 
(l2013h used we dges to measure th e RSD signal and reported / = 
0.719 zb 0.094. iReid et all (l2014l) has done the analysis at small 
scale 0.8 h“^Mpc < s < 32 h ^Mpc with halo occupation 
distribution model and Planck be st fit cosmology a nd measured 
fag{z = 0.57) = 0.45 zb 0.011. Ip^id et al.l (l2014h provides the 
strongest constraint on the growth rate but this analysis has signifi - 
cant modeling and cosmological assumptions. iMore et al] (l2014h 
measured the constraint on Dm and ag using a combination of 
abundance, clustering and galaxy-galaxy lensing. They have re¬ 
ported a constraint on fag by assuming the General Relativity lin¬ 
ear theory prediction for growth rate (/ = Our measure¬ 

ment is competitive with all RSD measurements from large scale. 


7 DISCUSSION 

We have presented an analysis of Redshift Space Distortion (RSD) 
using the SDSS-III BOSS DRll CMASS sample, and have mea¬ 
sured the monopole and quadrupole moments of galaxy auto cor¬ 
relation function at effective redshift of 0.57. We have used CLPT- 
GSRSD to model the Legendre moments of redshift space galaxy 
auto correlation function. The model used here does not work 
at small scales due to non-linearity, and measurement of corre¬ 
lation function from data shows systematic error at large scales. 
Therefore, we have adopted a conservative fitting scale between 30 
h“^Mpc and 126 h“^Mpc, which we chose with the aid of a suite 
of perturbation theory mocks. Our measurements of linear growth 
rate (fag), angular diameter distance (Da) and Hubble constant 
(H) at effective redshift of 0.57 don’t assume ACDM-GR evolu¬ 
tion by virtue of using Alcock-Paczynski parameters (ay, a^) in¬ 
dependent of cosmology at current epoch (z = 0). This approach 
makes these measurements suitable to test the predictions of vari¬ 
ous alternate models of gravity and c osmology. _ 

Our results are consistent with ISamushia et al.l (l2013h . who 
performed a similar analysis on the same data set. However, the 
perturbation theory models used in the two analyses are different. 
Our model (CLPT) performed better on N-body simulation com- 
pared to the Lagrangia n Perturbation Theory (LPT) model used in 
ISamushia et al.l (l2013h . . We have seen marginal improvement in 
the measurement uncertainty compared to the previous analyses. 
This is the first use of CLPT-GSRSD to measure both cosmology 
and growth s imultaneously from the galaxy redshift survey. It has 
been used bv iHowlett et al.l (l2015h to measure the growth rate with 
fixed cosmology for SDSS main galaxy sample. We couldn’t use 
our model at smaller scales because our mocks cannot be trusted 
in this range. In the future we may be able to extend this model to 
scales as low as 20 h“^Mpc if a reliable technique to test them on 
realistic mocks can be developed. 

The linear growth factor has been measured in many redshift 
surveys between redshift of 0 and 1. Our measurement provides an 
important data point to study the evolution of the linear growth fac¬ 
tor with redshift. The absolute value of fag and its evolution with 
redshift is quite sensitive to the model of gravity. These measure¬ 
ments will provide a good test of the general theory of relativity 
and the standard model of cosmology on the largest distance and 
time scales. It is possible to use these measurements to constrain 
flatness of the Universe and the dark energy equation of state pa¬ 
rameter. These measurements also have the ability to constrain the 
parameters of alternate theories of gravity and dark energy. 

The next-generation surveys are going to be even more power¬ 
ful, which will provide better measurement of correlation function 
and measurement of fag, hence better understanding of cosmol¬ 
ogy and gravity. The error in the measurement of the correlation 
function is much smaller at small scales, which has not yet been ex¬ 
plored in this paper due to our inability to test the theoretical model 
in this range. We can tap into the potential of small-scale clustering 
using RSD measurement when we can model the nonlinear cluster¬ 
ing at small scale either analytically or using fast simulations. 
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Osumi for providing systematic weighted correlation functions as 
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astro-ph preprint archive at arXiv.org. The analysis made 
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APPENDIX A: IMPROVING THE EFFICIENCY OF 
MCMC 

The most computationally expensive part of likelihood analysis is 
the calculation of CLPT correlation function and velocity statistics. 
In order to make this high dimensional optimization problem rea¬ 
sonably efficient, we have adopted two modifications. First we op¬ 
timized the sampling of k in the linear power spectrum used as the 
input to the CLPT perturbation theory. Second, we have discretized 
a small subspace of the parameter required for the perturbation the¬ 
ory in order to avoid doing almost the same calculation thousands 
of times. 


A1 Power Spectrum Sampling 

The CLPT’s runtime depends on the number of k points sampled 
in linear power spectrum. However, if we reduce the sampling of 
power spectrum too much, the integrals involved in the CLPT the¬ 
ory might not converge. Therefore, we need to minimize the sam¬ 
pling in k to reduce the calculation time but keep the sampling suf¬ 
ficiently high to avoid the convergence problem. We have run an 
optimization of the sampling in k by checking the convergence of 
the correlation function produced by CLPT and achieved the best 
case runtime of about 1 min for CLPT correlation function, which 
doesn’t have any convergence issues in the scale of interest, which 
is up to 130 h“^Mpc. 

We start with the finely-sampled power spectrum and kmax > 
100 Mpc“^h. The initial power spectrum (Figure |AT] shown in red) 
has linear sampling fork < 0.5 with 1030 points. The sampling is 
logarithmic for k > 0.5 with 500 points until kmax = 100. The 
CLPT takes 20 minutes to run with this sampling of power spec¬ 
trum. The correlation function converges until r = 200 h“^Mpc 
with the above mentioned sampling. We start with this power spec¬ 
trum and attempt to resample it in such a way that the run time 
decreases without introducing any error within the scale of inter¬ 
est which is assumed to be r < 130 h“^Mpc. The results of this 
optimization are summarized in Figure lAll 

Figure lAll shows the power spectra, which are vertically 
shifted for clarity. As we reduce the sampling of linear matter 
power spectrum, the computation time (shown in the legend of 
correlation function panel) decreases. As we reduce the sampling 
in k the correlation function fails to converge for large scale. The 
fastest runtime (Imin 18 seconds) produces a correlation function 
that doesn’t converge above 130 h“^Mpc. The corresponding sam¬ 
pling of linear matter power spectrum is as follows (shown in pur¬ 
ple in Figure I aTT) : 

k < 0.5 : linear sampling, 150 points 
0.5 < k < 1.8 : logarithmic sampling, 90 points 
1.8 < k < 10 : logarithmic sampling, 90 points 

A2 Discretizing a subspace of the full parameter space 

We wish to perform the likelihood minimization on the full param¬ 
eter space of Planck and RSD parameters. This approach creates a 
nine-dimensional parameter space, which requires millions of like¬ 
lihood evaluations; hence, the perturbation theory best case runtime 
of 1 minute is still too long to achieve convergence in a reasonable 
time with feasible computing resources.The CLPT, however, de¬ 
pends only on the input linear power spectrum which is the function 




Figure Al. Input power spectrum and output correlation function for four 
different run times . The power spectrum and correlation functions are verti¬ 
cally shifted in order to show the changes clearly. It is clear from the curves 
that as we sample the power spectrum more sparsely the correlation func¬ 
tion does not converge at large scale. This result decides our maximum fit¬ 
ting scale as 126 h“^Mpc for this analysis. The lower panel shows the 
fractional deviation of correlation function from the completely converged 
(red line) correlation function. The solid lines represents monopole and the 
dashed lines quadrupole. The black line demonstrates the fractional error in 
our measurement. The error in theoretical model is much smaller than the 
measurement error. 


of four cosmological parameters ris, Hq}. There¬ 

fore, we can significantly reduce the number of evaluations of the 
CLPT correlation function if it is done on the subspace of the full 
nine-dimensional parameter space. In order to avoid evaluating the 
CLPT repeatedly for the same cosmology, we have discretized the 
four-dimensional subspace and run the CLPT calculation only once 
for each grid in this subspace. Any repetitive call of CLPT for the 
same grid point will use the stored result from the previous eval- 
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uation for that grid point, Significantly increasing the speed of the 
entire optimization problem. 
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